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An approximate treatment of modes in a waveguide of variable width is presented. 
It is assumed that the boundaries satisfy impedance-type boundary conditions. The model 
consists of two parallel plate waveguide regions connected by a linearly tapered section. 
The results have application to the theory of VLF radio wave propagation when the iono- 
spheric heights are not constant along the path. 

1. Introduction 

The waveguide mode theory has been quite successful in explaining inau}^ of the broad 
features of propagation of VLF radio waves to great distances. Furthermore, the diurnal 
variation of phase of a cw carrier can often be interpreted in terms of the relative amount of 
solar illumination on the gi'eat circle path. In carrying out such an interpretation, it is usually 
assumed that a waveguide mode passes smoothly through the day/night transition without 
change of field pattern across the guide. Indeed, the experimental data show little evidence 
of distortion of the field in this transition region. Exceptions are indicated in the data of 
Lauter and Sprenger [1952] who observed amplitude variations in East Germany from the 
GBR transmission from England whicli indicated an undulating field. He explained this 
qualitatively in terms of diffraction from a ledge of ionization which w^as truncated at sunrise 
and sunset. 

An alternative approach is followed here which w^ould appear to be more appropriate for 
propagation to gi^eat distances. The model chosen is a parallel plate waveguide of constant 
width h which is connected to a waveguide of constant width h-{-Ah by a linearl}^ tapered 
section. The nature of the waveguide modes in the two parallel plate waveguides is well 
known. The character of the modes, in the tapered section, is first obtained by examining 
the permissible solutions in idealized wedge regions. Then an approximate metliod is applied 
to evaluate the coupling between the three waveguide regions. 

2. Wedge Solutions 

We consider a wedge consisting of two nonparallel plane surfaces. With respect to a 
cylindrical coordinate system (p,</),3), the apex of the wedge is the z axis and its surfaces are 
at = and (t>='^. The situation is shown in figure 1. For simplicity the fields are taken 




Figure 1. The wedge region. 
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not to vary in the z direction so that d/d2=0. This restricts the direction of propagation to 
be normal to the apex of the wedge region. Furthermore, the magnetic field is taken to have 
only a z component. Since the lower surface (0 = 0) is to represent the ground, it will be 
assumed that Ep = for 0=0 and 0<Cp<C°°- The upper boundary which is to represent the 
lower edge of the ionosphere is assumed to be characterized by a surface impedance Z. Thus 

Ep= -ZH, for 0=iA and 0<p< oo . (1) 

A suitable solution within the region 0<^(t)<^\f/ is of the form [Wait, 1959] 

H,=AZ,{kp) cos vet) (2) 

where A is a constant and Z^ is a cylindrical Bessel function of (complex) order v and argument 
kp. The choice of the factor cos v(j> automatically assures that the boundary condition at 
0=0 is satisfied. 

The boundary condition at ^^i/' can be rewritten in the form ^ 

^=-uo,pZ{p)H, (3) 

where it is indicated that the surface impedance is a function of p. To facilitate the solution, it 
is necessary to choose the p dependence such that 

Z(p)=Z,^ (4) 

P 

where po and Zq are constants. Thus 

~^=--ieo)poZoH, at (t)=\l/ (5) 

where e is the dielectric constant in the wedge region. On applying eq (5) to eq (2) we see that 

V tan v\[/^iecjopQZo. (6) 

Solution of this equation yields a set of values of v which are denoted Vg. 

We shall now demonstrate that these modes are orthogonal over the angular domain of 
the wedge. For example, consider 

Is,Q= I cos {vs4>) cos {v^<t))d(t) (7) 

where Vg and Vq both satisfy eq (6). Since 

J. _ sin (vs-v,)yp I sin (vs+Vg)^|^ .^. 

it is not difficult to show that 

which demonstrates the orthogonality. The normalization constant for the modes is 

It is now a simple matter to express the field Hz{p,^) at any point in the wedge in terms of 
its value over the surface p=pi, O<\0<CiA- Since the modes form a complete set, 

' The time factor is exp (icot) . 
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H^{p,ct>)^J2AZ.^{kp) cos vs<t> (10) 

where the summation is over all modes. Then, because of orthogoruility, the coefficient As 
for the sih mode may be expressed as 

p 

(11) 



where 



^ ^r sin (2.,^) -] 

Ps= Hz{pi,<t>) cos Z^,</>(/</). 



Thus, if we can find a suitable approximation for H.(pi,(f)) the corresponding function H^{p,(t)) 
can be obtained by performing an integration. 

Before proceeding further it is desirable to discuss something about the properties of the 
modes in the wedge region. Until the source itself is prescribed, the radial or p dependence of 
the field is not known. However, in general, the cylindrical Bessel function Z^sikp) may be 
represented as the linear combination of two IIank(^I functions, of [\\c first and scM'ond kind, 
of order v^. Thus a mode is of the form 



{ 






If kp >>1 and kp y>\vs\ it is permissible to replace the Ilankc^l functions hy tlieii' Dc^hye or 
second-order approximations [Wait, 1959]. Thus 

/:/lf (kp)^(-r-^. — Y (^xp [±ikp{sin a-a cos a)]^^^'^/^ (13) 

^ \Tkp sin aj ^ 

where a is defined by 

Vs=kp cos a. 
An alternative form is 

m' (kp) ^( ^^psi„ J ' exp [ ± ikp sin a] exp [ T i («''»+0]- (14) 

It is evident that these have the character of traveling waves propagating in negative and 
positive p directions, respectively. In this case the mode equation Vg tan Vs^=i^^PoZo can 
be written in the form 

Ctan kCh^ = iZ/r]Q 
or 



C+Z/rjo 
where 



(15) 



C= COS a=Vs/kp, 
and 



p 



h'=plk. 

It is seen that eq (15) has [)r(H*is(^ly the form of the modal equation for a parallel place wave- 
guide of equivalent height // and surface impedance (on one wall) equal to Z [Wait, I960]. 
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The choice of the kind of Hankel function (i.e., first or second) is dictated by the nature of 
the problem. In general, a linear combination is required. However, if the waves are excited 
over the aperture plane p=pi then for p>pi in an infinite wedge region the Hankel function of 
the second kind only is needed. 

3. Composite Problem 

We shall now consider the composite waveguide problem mentioned in the introduction. 
The situation is illustrated in figure 2. 

A Cartesian coordinate system, {x,y,z), is chosen such that the earth's surface (assumed 
flat) is the plane 2=0. An ionospheric reflecting layer of surface impedance Zi is located at 
z=h and x<C^O. A similar reflecting layer of surface impedance Z2 is located at z=h-\-Ah for 
x^b. The transition section extends from x=0 to x=b with a linear variation of height 
from A to h+Ah. The surface impedance of the surface of the transition section must be a 
linear function of p in order to make use of the orthogonal set of modes discussed above. It is 
convenient to choose Z=ZoPo/p where po=a sec x// which is equivalent to 



Z=Zq — , — for the interval 0<x<6. 

x+a 

A single waveguide mode is now considered to be incident from the left. It is taken to 
be of the form [Wait, 1960] 



Hz=e~^^''^v cos kCpZ 
where Sj,= {l — Cl)^ and Cj, is a solution of 



where :2^=0, 1,2,.. . and 






(16) 



(17) 



Ai= 



^Zi 

^0 



The integer jp is the order of the waveguide mode under consideration. A general incident 
field is obtained by superimposing all possible p modes. 

In a rigorous treatment of this problem we would need to set up a general representation 
of the fields in the three waveguide regions. To evaluate the unknown coefficients it is then 
required to match the tangential electric and magnetic fields across the aperture planes at 
x=Q and x=h. This leads to a very cumbersome (infinite) set of equations which must be 
solved for the infinite sets of coefficients. The problem is also complicated since the plane 
a:= constant is not a coordinate surface in the natural cylindrical coordinates for the wedge 
region. 

In order to obtain a useful result, albeit approximate, a number of simplifying assumptions 
are made. The justification for these are mainly on physical grounds. Essentially, the idea 
is to neglect the reflected waves at the aperture planes at x=0 and x=h. This approach has 
been used by Solymar [1959] who treated a number of analogous problems in microwave guides, 
and by Furutsu [1957] who considered related problems in radio propagation over mixed paths. 




Figure 2. The transition 'problem {not drawn to 
scale since \p^l). 



x = o 
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In view of the above assumption, the field in the wedge transition region has the form 

H,=^AMf(fcp) cos v,ct> (18) 

where Vs is a sohition of 

Vs tan Vs\p=ieo:poZo. 
The coefficients As are given by 

where 

P,= I H,{po,(l)) cos VsCfdv, 

If the length of the transition section, b, is large compared with the heiglit change, Ah, it 
is reasonable to assume that H^{p(^,4>) is given adequately by tlie incident mode. This has the 
form 

j^inc^ ^-ikxs^ COS kCj,z (20) 

and thus 

H, (po,0) ^e-^'^^'v cos (kpoO,,ct>) (21) 

where 

2^o=(Po— 2^)i — Po COS xp 



Therefore 



^^''-£ cos (kpoC,<t>) COS (.,</>) exp [-^fcS,po (f ^f )] ^^</> (22) 



where the superscript p is to indicate that this is the value appropriate to the 7>>th incident 
mode. The function P^f in the form given above niay be rewi'itten in the form 

p(.) _ e-^kSvf>o^'imkpoC,,Vs,kpoSJ2,rP] (23) 

where the integral / of the four variables, is expressible in terms of Fresnel integrals (see 
appendix). This is restricted only by the previously imposed condition that \^<^<^1. 
A somewhat simpler approach is possible if, in addition, kpo\l/^<^<Cl. Then 

pjP)^e-iA-poS^^2/2 / p (.^g ^^ COS v<t>d<t>+^^^^^^ r </>2 COS M</> COS pct>dcl>\ (24) 

where iJi=kpoCj,^khCp/\p and v=Vs' In view of 

V tan vyp'^ieciopoZQ 
and 

jjL tan /xi/'=?*€a;poZi 
it follows that (when Zq=Zi) 



cos ti4> cos v4>d(i>= -< 2 L 2^^ J 



-7— I tor i'=ju 



(25) 



for vy^p.. 

The second integral of eq (24) can be reduced in the following way 
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</)2 COS /x0 COS v(l)d<l)^^ 0^ COS {^i—v)(t)d(t)+:^ </)^ cos {ix + v)<l>d(t> 

=?r7 ^ •'^^ cos X^X-|-;r7 ^ — r^ X^ COS xc?x* 

2(m-^)^ Jo 2(/x+v)3 Jo 



^ [2(/x-^)iA COS {^-v)xP+[{fx-pyxl^'-2] sin (iu-^)\^] 



1 



[2{^+v)xP cos (m + .>+[(m+^)V'-2] sin (m+^¥]. (26) 



' 2(m+v)3 
This latter result is usable for /jl^v. When fx = v the integral may be handled in the manner 

(t>^ cos v(t)d<t)=7z 0^(l+cos 2jLi(^)rf</)=— +7^^-^ x^-cos2xdx 

^ Jo o (Zfi) Jo 



■ 6 "^8m' 



=^+^3 [(W) cos (2m^) +[2(m^)^- 1 ] sin (2m^)J. (27) 



While this is essentiullj' the sohition of the problem it is desirable to make some fm'ther 
simplifications to permit discussion. Since the important modes at VLF are near grazing 
we ca]i approximate m and v as follows [Wait, 1960] 






Thus 

r^ - for p=s 

I cos ii(t) COS v(t)d4>= < ^ 

^' lo for p9^s. 

I 6^ cos U0 cos ^</)C/0 = 7 r^ COS iriv — Sj + 'z , — r^ COS iriv^s) 

Jo (m— ^) (m + ^) 

On the other hand, if ^=^, the integral is given by 

J^ 0H.os.0:/^4-^-,^^, (29) 

We are now in a position to write down some relatively simple formulae which describe the 
transition. First, we see that 

pr=\[^+o{r)] (30) 

where 0(^^) is a quantity which is of second-order smallness in yp. Consequently 
In the case of the modes where p9^s, 
and therefore 

„_iW!ei)!:!r_j__ I 1 ,33. 
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A measure of the mode conversion is expressed by the ratio 

Wlien the incident mode is of order p=l, tliis ratio is simply expressed by 

Pl^'^imr i-2s -1 

P.<"= x^ Ls'(s-l)'J 
and for s=2 

IT _ikH 3 

We are now in a position to calculate the field in the parallel plate region for aC> h. Here 
the field may also be written as a sum of modes in the manner 

^z=Z;j5»e-'*'^-»^"cos ()t<7„3) (35) 

n 

for 0<^.c<A+AA. The modal equation in this ease is 

C;tan[/:r,(/^+A/^)] = iZ,/r;o 
and in view of tlie orthogonality of th(» modes 

lH,{x=b)] cos (kCnZ)(h 
B„=^^^ ^ ^ , (36) 



, , sin 2kCJi 

1+ 1^ 



2kCnh 

where n = A+A^. Following earlier reasoning, the field H^ over th(^ aperture x = h is now 
approximated by the field incident from the wedge transition region. Thus 

H,{x=h)^H,{p',,<i>)e'-''^''^ (37) 

where 



Po=[h''+{h+m'\ 



1/2 



4=(po—^^)^— Pottos ^ 
'\2 2) 






and 

H,{p',,<i>)^^Al^^mf{kp',) cos v,<t>. 

s 

Again, for the sake of simplicity, one may choose Z{b) = Z2 so there is no discontinuity in tlie 
surface impedance. 

4. Summary of Final Results 

We are now in a position to summarize and restate the results obtained above. The fields 
in the three waveguide regions have the following (approximate) forms 

H,=e-'^^> cos {kC'^z) for x<0, 

ffz^( T , , j ^ ocs^ 'Jo ^""^ co^[kCs{x)z] forO<a:<6, 

H,^(^j^y S e,e~''''^' ^'~'^ cos (kC:;z) for x>b. (38) 
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The respective modal equations are 

C tan {kC'h)=iZilr}o, 

C{x) tan [kC{x)h(x)]=iZ(x)/vo, 

C tan kC''(h+Ah)=iZ2/vo- (39) 

The coefBcients ag and ^^ are expressible in terms of Fresnel integrals which are functions of the 
angle ^, b, Zi, Zo, and Z{x) as discussed in previous sections. Furthermore, in order to assure 
orthogonality of the modes in the wedge region, it was necessary to choose 



Z{x)=Zo 



a+x 
and 

h(x) = h+TAh 



>- for 0<a;<6. (40) 



It can be expected that the nature of the coefficients a, and 13 g are not significantly influenced by 
the finite values of the surface impedances Zi, Z{x), and Z2. This is particularly so in the case 
of grazing modes where the phase shift on reflection at the upper boundary is near tt. Therefore 

and 



^^^^^^^[l^-ip+l+ir] (-1)'"^ 



for 
If 

and 






for 

5=2,3,4. . .. (41) 

Now at the junction at x = b, the situation is very similar. Thus I3p^l and fig^—as for 
8 9^ p. This implies, of course, that reconversion from modes of order s{7^p) back to the mode 
of order p are neglected. In view of our previous assumption about the smallness of y//, this 
is entirely justified. 

We see from the present analysis that the mode conversion is probably an extremely small 
effect and can be neglected in most cases of practical interest. Indeed, experimental data on 
amplitude and phase variations do not appear to have any pronounced and reproducible fea- 
tures which could be attributed to mode conversion at a day/night transition.^ However, 
local ionospheric disturbances may well produce abrupt height changes which could lead to 
effective values of \f/ which are not small. Further effort is being directed toward a better 
understanding of these phenomena. 



I thank Mrs. Alyce M. Conda for pointing out a number of copying errors. I am also 
indebted to Prof. J. E. Lindsay for a number of valuable comments. 



2 A. Q. Jean (personal communication) indicates that there is a consistent phase distortion for the path NBA (Panama) to Boulder. It is pos- 
sible this could be attributed to the interaction with the spurious second mode produced at the transition boundary- 
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5. Appendix 
Evaluation of 

I {fi,v,a,\l/)= COS A10COS v(j)e^'^'*'^d^. 



where 



The above integral can be transformed in the following manner 

1 r^ ■> 1 r^ 

- cos (ju— J')</>6*«*^^6;<^+^ COS (M + »')e'«*^V<^ 

4 Jo 4 Jo 

4 Jo 

^e '" \ e ^ '^^ d<l>+ . .. 
4 Jo 

-(t)X*+'£^)- 



4 

_ /? 4a 

4^ 



and 



Thus 



.,..„.)=qK {'i©'(*+'£')]+'i(v)'(*-^-')]} (0 



-4^ {"{(l)'(*+'^)]+'{(l)'(*-'^)]} ©' <«' 



where F{x)= I e ^ 6?a^ is a Fresnel integral. 
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